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WAVE GROWTH IN A STRONGLY TURBULENT PLASMA 

ABSTRACT 

We derive an equation for the average non-linear growth (damping) rate 
y of an electrostatic mode with wave vector k for an ensemble of electro- 
statically turbulent, unmagnetized plasmas. Strong turbulence alters particle 
orbits during the course of wave growth and hence brings particles in and out of 
resonance with the wave. The wave growth rate in a strongly turbulent plasma 
therefore depends on the initial velocity distribution of particles not only at 
u = k • v = a k /k , the wave speed, but in the vicinity of this speed as well. We 
consider two orbital effects: turbulence modification of ensemble average orbits 
and turbulence produced orbital dispersion about the average orbits. The first 
of these leads to a shift in the central u of the resonance. Both contribute to a 
resonance broadening. 
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WAVE GROWTH IN A STRONGLY TURBULENT PLASMA 


INTRODUCTION 

The diffusion of resonant particles in the presence of strong plasma turbu- 
lence was the subject of a previous paper 1 (hereafter referred to as I). We con- 
sidered in I an ensemble of plasmas and found that strongly turbulent fluctuations 
modify the statistical orbits of resonant particles and hence affect the ensemble 
average diffusion rate. The fluctuations (1) modify the ensemble average orbits 
of resonant particles and (2) introduce dispersion in orbits about the averages. 
Both effects are quadratic in the magnitude of the fluctuations. 

In strong turbulence the coherence time between resonant particles and 
turbulence depends on the strength of these two effects. Initially resonant 
particles may be driven from resonance, or conversely initially non-resonant 
particles may be made more resonant by the turbulence. The more rapidly that 
these processes occur, the more significantly is the diffusion tensor altered. 

Turning now to the topic of wave growth (damping), we recall that in the 
Landau 2 analysis y k , the growth (damping) factor for an electrostatic wave of 
frequency and wave number k in an unmagnetized, one dimensional electron 
plasma is 
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( 1 ) 



proportional to the slope of the electron one particle distribution function f at the 
resonant velocity v = a k /k (oj s is the electron plasma frequency). Tacit in 
Equation (1) is the assumption that resonant electrons move in straight lines 
over the time interval !r u i ' 1 . Large amplitude fluctuations, as may occur in an 
unstable plasma, lead to orbit modifications which violate this straight line 
approximation. 

In the strongly turbulent regime we consider an ensemble of three dimen- 
sional plasmas just as in I. We seek here the value of the ensemble average 
growth rate <(> k ^ for an electrostatic wave of frequency a k and wave vector k. 

Both statistical orbit modifications discussed in I also affect the form of 
k )> . Because the average velocity of each particle is changing and its velocity 
is also dispersing during the course of wave growth, a wave of phase velocity 
a R /k interacts not just with particles whose parallel velocity u = k • v equals 
a k /k at the beginning of growth; it interacts also with particles whose initial u’s 
lie near a k /k and which thus have a high probability of resonating (i.e., having 
u = a k /k) during the growth period . 

In a strongly turbulent three dimensional plasma thus depends on the 

value of <f !>over a band of u. Here <Cf ^is the ensemble averaged one particle 
distribution function integrated over velocity components perpendicular to k. If 
the wave-particle interaction were the same for all resonant particles, one would expect 
<(y k )> to be of the form 
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the distortion of the resonance being of second order in the fluctuation ampli- 
tude SE and having contributions arising from both average orbit modification 
and orbit dispersion. 


In this paper we calculate <y k > for electrostatic waves in a spatially homo- 
geneous, three dimensional, multi-component plasma. The turbulence which 
deforms particle orbits is likewise electrostatic. (Indeed, the mode whose 
stability we are considering is one of the waves comprising the turbulence.) The 
importance of average orbit distortion and orbit dispersion here lies in the fact 
that these effects can act as an amplitude-dependent stabilizing factor for waves 
which are linearly unstable 


[ <?,> (ir) > 0 for T " 0 

Dupree 3 has already included the effect of orbit dispersion in his one 
dimensional calculation. His pioneering work has been amplified and applied to 
specific linearly unstable modes by several authors . 4 ' 12 We here add the effect 
on the growth rate of average orbit alteration, an effect, like dispersion, quad- 
ratic in SE. 


As in I our approach is via Dupree's 3 statistical orbits technique. We 
make frequent use of results derived and approximations discussed in I. The 
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reader who is interested in the details of this calculation is thus advised to 
study the papers in conjunction with one another. 

CALCULATION OF (r h > 

In each realization of the ensemble, Poisson's Equation relates the fluctu- 
ating electrostatic field 8 E to the fluctuating charge density 8 p and thus to § f . , 
the fluctuating part of the one particle distribution function for each species i 


V‘SE=47rSp = 4 


v 




d 3 v S f (x, v, t) 


( 2 ) 


In the weak coupling approximation to strong plasma turbulence, 8 L is in 
turn related to its initial value 8 f. (x, v, t 0 ), the ensemble average one particle 
distribution function <f>. (v, t) for species i, and 8E(x, t) itself by 


8 f. (x, v, t) = <U. (t, t Q )> 8 f. (x, v, t Q ) 


-- r 


( 3 ) 


B <f>. 

d r <U. (t, r)> 8 E (x, r) • — — — - (v, r) 

O V 


The propagator ^U^t, t 0 )^> is the ensemble average Vlasov propagator for 
particles of species i. <Cu(t,t o ))> is defined by Equations (9) and (10) of I and 
contains the effects of the turbulent electrostatic waves on particle trajectories. 
The relationship between <U(t,t 0 ))> and the characteristic phase space trajec- 
tories of plasma particles is developed by Birmingham and Bornatici. 13 The weak 
coupling approximation is discussed in detail in I. 
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We next plug Equation (3) into Equation (2). The initial value contribution 
serves as a driving source for the waves. Since we are here concerned only with 
the dielectric properties of the plasma, we assume that S f^x, v, t Q ) is a well 
behaved function and simply represent its effect in Poisson's Equation by 
S(x, t; t Q ). Thus 


V- SE = 



J dr <U.(t, t)> 8E (x, r)* 

t n 


XOi (V, T) 


+ S (x, t; t.) ( 4 ) 


where 

4 7T n. q? 

= 1__1 

1 m. 

1 

is the plasma frequency associated with the i th species. 

In Equation (4) •\U i (t, r )^> operates on both SE and B<f > i /Bv . We shall, 
however, assume (1) that the scale length of <f >. in v is large compared to the 
incremental change in velocity coordinates of particles in the growth time 
|-y k | 1 of any mode and (2) that the time scale on which <U>. changes is like- 
wise long compared with |y k | _1 . [In the important resonant particle region of 
phase space, < f> 4 evolves according to the diffusion equation derived in I on the 
time scale of several correlation times r c .] Corrections introduced by includ- 
ing the finite velocity and time scales of <f_> A are of the same order of magni- 
tude as those encountered in considering higher order terms in the weak coupling 
approximation. 
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A cumulant expansion of <U.(t, r))> 14 - 15 (cf. Equation 14 of I) puts Equation 
(4) in the form 

V- SE = -2> ? j V f exp j^ x i C 7 ")/* ‘^+ [<A X* (r) A x* (r)^> 


- <f A x* ( T )> <C A X i ( T )> 


1 B<f>. 

: V Vj 8E (x, t) * — ^ 1 + S (x, t; t Q ) 


x + < A x* (t) ) is the ensemble average phase space position coordinate at time 
r of a particle of species i which in each realization of the ensemble is at the 
point x, v at time t. Similarly, 


I = [<A x* (r) A x* (t)> - <A x* (r)> <Ax* (r)>] 


is the position dispersion tensor for such a particle. 

We have derived in I expressions relating <^A x* (r))> and I to the particle 
diffusion tensor D for an unmagnetized, spatially homogeneous, and quasi- 
stationary (in the sense that |y k | -1 is longer than the cut-off time 
of autocorrelation functions), but strongly turbulent plasma. For particles of 
species i these expressions are 


<A X* (T)> = - V (t - r) - {1 ~ T) JL Tr (g. (V, t - r)} 


( 6 ) 


h - T (t ~ T ^ 3 = i ( V> t - T ) 


(?) 
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In Equation (6) Tr {D.(v, t - r)} is the trace of the matrix of D.. The diffusion 
tensor is essentially zero for all particles except those resonant with the turbu- 
lence. For resonant particles D i is time dependent, building up from zero until 
the time when fluctuation correlations (as observed by the resonant particles) 
have died out. Beyond this time D i has the asymptotic value given by Equation 
(28) of I. 


Let us now introduce the forms Equations (6) and (7) into Equation (5) and at 
the same time Fourier transform in x: 


(t) 
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Tr {D. (v, t - t)} - 1 (t - r) 3 k k: D (v, t - r)^- S E k (r) 

KOi i S, (t ; t 0 ) 


B v 


k 2 


(8) 


Asymptotically in time (t -» -oo) we expect only the natural oscillation modes 
of the system to persist and write SE k (r) a exp -i <(« k ^> r. The ensemble aver- 
age eigenfrequency <(a) k ^> occurs because in Equation (8) the dynamic response 
of the plasma to all waves has been averaged over and represented by the Dj ’s 
We thus obtain the dispersion relation 


/» /* 00 

1 =_ l - jyi I d 3 V d z exp -fi T( <a) R > - k-v) z 
k i 0 

Tr <g i (v, ,»]. - | z 3 k k: D. (v, z)}k' 


(9) 
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For slowly growing (decaying) waves, 


I <r k > I = |lm«o) k >)| « | <a h > | = |Re(<<^»|, 

we obtain by taking real and imaginary parts of Equation (9): 

1 = - -L y ' a? Id 3 v k • — — — 1_ I d z s in - k • v) z - z 2 A £ ( z)] exp - z 3 C (z) 

(10a) 


Oi,> 


r— , f 3<f>. f* 
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d 3 v k * — - — | d 2 cos [(<^ a k > - k • v) z - z 2 A. (z)] exp - z 3 C A (z) 


^ f 3 < f ) . f 00 

^ / w? d 3 vk 1 ~ d z sin [( < a-> - k ■ v) z - z 2 A. (z)] exp - z 3 C { (z) 

B < a k> “ 1 J " a v 1 


(10b) 


We have introduced here the definitions 


Ai ( z ) = ik*-_L Tr {D x (v, z)> 

2 d v « l 

and 

Cj O) = ik k: D. (v, z) 

for notational convenience. Equation (10a) determines </y k ^> which when inserted 
into Equation (10b) yields Kp-^> • 

We assert that in Equations (10) A { and C A may be neglected in the integrals 
involving the sine function. The main contributions to these integrals comes 
from the non- resonant region of velocity space which is not greatly affected by 
the turbulence. Equation (10a) thus reduces to 
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a k - k v 


the form obtained in linear theory. (Since a k is not affected by the turbulence 
we drop the < > .) 


On the other hand resonant particles play the dominant role in the integral 
appearing in the numerator of Equation (10b). Recall that 


1 r 

im 

i 77 / 

Ho ♦'o 


d z cos [(a, - k • v) z] exp - d z = 1 : 


d-* 0 77 d 2 + (a - k- v) 2 


is one representation of 8 (a h - k • v ). We thus obtain 




d 3 v k 


3 <f>; 


d z cos [(a h - k • v) z - z 2 Aj (z)] exp - z 3 C £ (z) 


Ok> =- 


d 3 v k'3 (f)./3 v (a k - k • v)” 1 


Note that for an electron plasma with no streaming, Equation (11) predicts 
= ± co and that in the absence of turbulence Equation (12) then yields 


< 7 k > = 


77 a k 'e y, H 


(?) 


in agreement with Equation (1). 


As in I the strongly damped nature of the z-integrand in Equation (12) allows 
us to approximate A. (z) and C. (z) by their small z expansions 
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(14) 


A. (z) = 1 z 3 k-T" 1 k' (a k « - k' • v) ^|S E k » | 2 > = z 3 A. 

6 m 3 V 

C 4 (z) ^kk:7> k* <|8E k .|- 2 > 

m i k' 

(15) 

1 q? -v 

- - z — k k: <8 E (x, t) 8 E (x, t)> = z C. 

3 m 2 


Here 8 < is the amplitude of the electrostatic turbulence mode with wave 

vector k ' . 


We now plug the forms Equations (14) and (15) into Equation (12) and obtain 


L f 3 <f>. P 

co? I d 3 v k - — I d z cos [(a k - k ■ v) z - z 5 A.] exp - z 4 


<*> 


3/3 


«*£«? pJ d 3 v k • 3 <f >. 


(16) 


/3 v (a k - k • v) _1 


Equation 16 is the main result of this paper. Effects of the strong turbulence 
are explicitly manifest: alteration of ensemble averaged orbits is the source of 
the z 5 - term in the cosine function and dispersion in orbits about the ensemble 
average results in the z 4 - exponential damping. 


The non-linear wave-particle interaction for species i is contained in the 
function 


G 4 (k, 


CU 


= k • v, 





z cos 


[( a k 


k u) z - z s Aj ] exp - z 4 Cj(17) 
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which in the absence of turbulence is just 77 S(a k - ku). The effects of turbulence, 
as we shall see, are to broaden this resonance and to shift its peak from u = a k /k. 

Note from Eqs. 14 and 15 that A* has dimensions (time) - 5 and has 
dimensions (time) -4 . In order of magnitude ~O(k/k 0 r 4 . r, ) and c i ~0(k 2 /koT 4 .) . 
where k 0 is a typical wave number of the turbulent spectrum. r c has been de- 
fined in I as the autocorrelation time of the fluctuations as observed by a reso- 
nant particle moving as if there were no fluctuations. ( r c is approximately the 
same for all resonant particles, irrespective of species.) Also from I, t 2 is 
the time for a resonant particle to diffuse (because of the turbulence) the distance 
2 77/k o (t 2 i is t 2 for particles of species i.) 

As stressed in I the weak coupling approximation to strong plasma turbu- 
lence limits us to a consideration of the case r c < < r 2 . The ratio Ai 5/4 /C i is 
hence in order of magnitude greater than unity. We conclude therefore that 
average orbit distortion is an important factor in determining the nature of the 
resonance function G i; since over the decay time of the exponential the z 5 factor 
in the cosine function has produced a considerable oscillation. 

The fact that A^ 7 vC^ _ 1 prohibits us from expanding the integrand in 
Eq. (17) with respect to small A*. Consequently there is not much that can be 
done analytically with the resonance function G j. To see the qualitative effects 
of finite Xi and C- we have therefore performed a numerical integration of 
Eq. 17 for a one dimensional example. The turbulence spectrum that we have 
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chosen is the same one discussed in I except that we here have taken Ak/k 0 = l/2 
and e = 10“ 5 . 

Plotted in Figure I is G i in units of a k against v = kv/a k . Recall that a k 
and k are the frequency and wave number of the mode whose stability is being 
studied. In Curve I we have taken Aj. = 0 so this trace illustrates the resonance 
broadening about v = l due to orbit dispersion as predicted and discussed by 
Dupree. 3 In Curve II we have added the effect of average orbit distortion. The 
resonance is altered asymmetrically: its maximum is shifted to v > 1 and it 
acquires an additional broadening primarily on the large v side of maximum. 

The height of the maximum Gi is also diminished. Were we to replace the 
resonance function Curve II by a rectangular function (in the spirit of Dupree's 3 
treatment of all resonant particle interactions with the turbulence on an equal 
footing) the rectangular resonance function should thus be asymmetrically posi- 
tioned with respect to v = 1. This result is the basis of our assumptions in the 
Introduction of a shift in the "center of mass" of the resonance to 
v = (\A) + (gj - g 2 )/2. 

Recall now that a resonance function of the form Curve II for each species 
when folded with 3(f). /3v into a v integration determines the sign of < y k > 
and thus the stable or unstable nature of the plasma. Near marginal stability 
< 7 k > * 0 and hence we here expect that the deviations from asymmetry in the 
G/s are decisive in determining the stability character of the plasma. 
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SUMMARY 


We have derived an expression for < y h > , the average non-linear growth 
(damping) rate of an electrostatic mode of wave vector k , for an ensemble of 
electrostatically turbulent, unmagnetized plasmas. The non-linearity arises 
because in the course of wave growth turbulence alters the orbits of particles 
from straight lines. The non-linearity manifests itself as both a shift in the 
central u = k • v for the Landau resonance and a broadening of this resonance 
in u. 

Two characteristic times enter the problem: t., the correlation time of the 
fluctuations as observed by resonant particles moving as if there were no 
turbulence, and r 2i , the time for a resonant particle of species i to diffuse one 
typical wavelength 27 r/k 0 of the turbulence. In this theory we consider the weak 
coupling approximation to strong plasma turbulence so that r c << r 2i . Average 
orbit distortion then becomes a significant factor in determining the form of <r k > 

Average orbit distortion leads to both a shift in the Landau resonance and 
a broadening of its width. We feel that the resonance shift is especially im- 
portant when the modes are linearly near marginal stability. 
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Figure Caption 


Figure 1. 


The resonance function G i with orbit dispersion (Curve I) and both 
orbit dispersion and average orbit distortion (Curve II). 
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